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ABSTRACT 

Aims. We compute photometric redshifts in the fourth public release of the Canada-France-Hawaii Telescope Legacy Survey. This 
unique multi-colour catalogue comprises u* ,g' ,r' ,i' ,z' photometry in four deep fields of 1 deg^ each and 35 deg distributed over 
three Wide fields. 

Methods. Our photometric redshifts, computed using a template-fitting method, are calibrated with a large catalogue of 16,983 high- 
quality spectroscopic redshifts from several spectroscopic surveys. 

Results. From the comparison with spectroscopic redshifts, we find a photometric redshift dispersion of 0.028 and an outlier rate of 
3.5% in the Deep field at i'^g < 24. In the Wide fields, we find a dispersion of 0.036 and outlier rate of 2.8% at i'^^ < 22.5. Beyond 
i'j^g = 22.5 in the Wide fields the number of outliers rises from 5% to 10% at i'^^ < 23 and i'^g < 24 respectively. For the Wide 
sample, we find the systematic redshift bias keeps below 1% to < 22.5, whereas we find no significant bias in the Deep field. We 
investigated the effect of tile-to-tile photometric variations and demonstrate that the accuracy of our photometric redshifts is reduced 
by at most 21%. We separate stars from galaxies using both the size and colour information. Applying our star-galaxy classifier we 
reduce the contamination by stars in our catalogues from 50% to 8% at i'^g < 22.5 in our field with the highest stellar density while 
keeping a complete galaxy sample. Our CFHTLS T0004 photometric redshifts are distributed to the community. Our release include 
592, 891 {i'^g < 22.5) and 244, 701 (i'^g < 24) reliable galaxy photometric redshifts in the Wide and Deep fields, respectively. 

Key words. Photometric redshifts - cosmology - galaxies 



1. Introduction 

It is now evident that the exploration of large scale struc- 
ture and the high-redshift Universe with the Canada France 
Hawaii Telescope Legacy Survey (CFHTLS) requires precise 
magnitudes and redshifts for millions of sources ('Be njamin et al j 
2007, Parker et al. 2007, McCracken et al. 2008, Fu et al.ll2008l 
Berge et al.ll2008t ,Kilbinger et al. 2008 TerenoeLaL 2008 ). 



To date, only "photometric redshift" techniques can provide 
(with a comparatively modest expenditure of telescope time) 
redshifts of enormous numbers of galaxies with sufficient preci- 
sion (~ 1-5%) to the faintest limiting magnitudes of the CFHTLS 
cosmological surveys. The construction of well defined, accu- 
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rate and reliable photometric redshift catalogues is therefore an 
indispensable task following photometric catalogue production. 

Several photometric redshift codes are now publicly avail- 
able and have been applied with reasoi iable success to 
many photometric catalogues of galaxies (see lHildebrandt et al.l 
(l2008h . and references therein). Some recent p hotometric 
redshift studies, like the COMBO- 17 survey (IWolf et al.l 
20031 , CFHTLS (Ilbertetal. 2006, hereafte r 106), SWIRE 
(Rowan-Robinson et al. 120081) . or COSMOS (iMobasher et al.1 
2007, Ilbert et al. 2008) . contain up to 1 ,000,000 galaxies as faint 
as i ~ 25. In particular, the "Le Phare" pho tometric redshift 
code (lArnouts et all 1 19991 lArnouts etani2002l 106) has shown 
to be well adapted for joint photometric and spectroscopic sur- 
veys like the CFHTLS. 106 used Le Phare with the CFHTLS- 
Deep photometric catal ogues and VVDS spectroscopic redshifts 
(iLe Fevre et al.l l2005bh to calculate photometric redshifts with 
an accuracy of ~ 3% at / < 24. The redshift distribution of 
sources has been used to calibrate the absolute gravitational 
shear signal presented in lFu et afl (|2008|) . 

By combining photometric and spectroscopic measurements 
they were able to derive a set of optimised Spectral Energy 
Distributions (SEDs). By computing the mean difference be- 
tween magnitudes in each filter of objects with known redshifts 
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and those derived from these optimised SEDs one can "tune" the 
CFHTLS photometry. This method improves photometric red- 
shifts for all galaxies, even those for which there are no spectro- 
scopic measurements. 

More spectra and deeper data can lead to a better calibration. 
With the advent of new CFHTLS photometric data and much 
larger spectroscopic catalogues we can now extend the applica- 
tion of Le Phare to more complex surveys composed of many 
MegaCam fields like the CFHTLS-Wide. 

Most photometric redshift studies either explore shallow 
very wide fields covering thousand of deg^ (see Ball et al. 2008 
for the SDSS and references therein) where visible photometric 
data are sufficient to sample the whole redshift range of galax- 
ies, or those which focus on deep beams of few deg^ (like the 
CFHTLS Deep or COSMOS) comprising both visible and near 
infrared photometric data. To date, the only moderately deep vis- 
ible survey currently available for photometric redshift studies 
covering a large area is the CFHTLS Wide. 

The aim of this work is the calibration and production 
of a flux limited photometric redshift catalogues based on 
the CFHTLS Deep and Wide surveys. When completed, the 
CFHTLS Wide will cover 170 deg^ in five optical filters spread 
over four separate regions of the sky which also contain subsets 
of several deep spectroscopic surveys. Ultimately, a catalogue of 
more than ten million galaxies down to /^^ - 24.5, with reliable 
five-band photometry will be available for photometric redshift 
measurement. In this paper, we use Le Phare to compute photo- 
metric redshift for the "T0004" CFHTLS release. In addition to 
the four CFHTLS Deep fields that were already analysed in pre- 
vious works, T0004 includes a new large Wide catalogue, cov- 
ering 35 deg^ in five bands u*,g',r',i',z' in three independent 
fields. 

To construct our photometric redshift catalogues we follow 
the method described in 106, but with some important changes. 
To calibrate the photometric redshifts new spectroscopic dat a 
are added to the VVDS "deep" sample (iLe Fevre et al.ll2005bh . 
It includes the D EEP2 redshift survey ( Davis et al.l l2003l"and 
Davis et al.l l2007l) . the new VVDS "wide" sample (Garil U et al.1 
,2008) and the "zCOSMOS" sample (Lifly et al. 2007). We first 
compute the new Deep field photometric redshift are then com- 
puted first and compared to the previous 106 Deep study. Once 
validated, we extend the analysis further to Wide data. We then 
challenge the photometric redshift catalogues against the new 
spectroscopic redshift samples to assess the robustness of the 
calibration, to derive a detailed error budget for each photometric 
redshift and to estimate the sensitivity of the method to param- 
eters like redshift or magnitude. Finally, we derive the redshift 
distribution and its field-to-field variance (cosmic variance). 

The paper is organised as follows. In Section lZTI we describe 
the data used in this analysis. We show how the "T0004" pho- 
tometry is used and we explain how the spectroscopic redshifts 
are evaluated in the WDS F02 and F22 fields and the DEEP2 
survey. The next Section describes the principles of Le Phare, 
and discusses recent improvements for the photometric redshift 
caUbration using new spectroscopic redshift samples. Sections]?] 
and ]5] focus on the photometric redshifts of the Deep and Wide 
CFHTLS sample. We then discuss the photometric redshift accu- 
racy and their robustness and reliability. The redshift distribution 
and its variance are discussed in the last section. 

Throughout the paper, we use a flat lambda cosmology 
(f2,„ = 0.3, Qa = 0.7) and we define h = //o/lOO, with 
H(i - 70 km s ' Mpc Magnitudes are given in the AB sys- 
tem. Photometric and spectroscopic redshifts are denoted by Zp 
and Zs, respectively, and Az represents Zp - Zs- 



2. Data 

2.1. Photometric data 

2.1.1. The CFHTLS 

The CFHTLS, a joint Canadian and French survey, is composed 
of three imaging surveys of different depths, shape and sky cov- 
erage: the CFHTLS Deep, the CFHTLS Wide and the CFHTLS 
Very Wide. The survey started in 2003 and will last for 450 
nights until finishing in January 2009. When completed, the 
CFHTLS Deep will cover 4x1 deg^ at a limiting magnitude 
of i'^g = 27.5 (point source, S/N=5, 1.15" aperture, seeing 0.8") 
and the Wide will cover 170 deg^ over 4 fields with a limiting 
magnitude of i'^^ = 25.5. 

The CFHTLS is conducted at the 3.6m Canada-France- 
Hawaii Telescope, equipped with the new MegaPrime prime- 
focus assembly and the MegaCam wide field camera 
dBoulade et alj]2000 ). MegaCam is composed of 36 2080x4644 
pixel EEV CCDs. MegaCam has a pixel scale of 0.186"/pixel 
and therefore covers the entire 0.96x0.98 deg^ unvingnetted 
MegaPrime focal plane. The CCD assembly leaves small gaps 
between each detector and two blind lanes of 62" width along 
the X-axis (E-W direction). For this reason each field is observed 
in a dithering pattern to fill the gaps between the CCDs and help 
in removing systematic CCD features. 

MegaPrime is equipped of u*,g',r',i',z' broad band filters 
that provide continuous coverage over the whole spectral range 
3500A < A < 9400A. Only the CFHTLS Deep and Wide sur- 
veys which are observed in the five MegaCam u*,g',r',i',z' 
broad band filters are suitable for photometric redshift compu- 
tation. 

The positions of the four deep fields and the three wide fields 
used in this work are listed in tableU] Each Wide field has a dif- 
ferent geometry and sky coverage: Wl (8° x 9° "tiles" of 1° x 1° 
each), W3 (7° x 7°), W4 (a 7° x 3.5° strip along the South- 
East/North- West direction), and is composed of a contiguous 
mosaic of MegaCam fields. The identification name and posi- 
tion of all fields can be found on the Terapix web pagefl 

Wl and W3 overlap with Dl and D3, respectively, while W4 
located a few degrees apart D4. Each MegaCam "tile" (which 
is an elementary 1° x 1° MegaCam area paving the sky) com- 
posing a Wide field overlaps with its four neighbouring fields. 
The width of overlap regions is around 2' in both directions. The 
overlap between pointings can be used for the astrometric and 
photometric calibrations. 

2.1.2. The T0004 release 

The data used in this work are part of the CFHTLS T0004 release 
produced at the Terapix data center. The release consists of all 
CFHTLS Deep and Wide images observed from June 1st, 2003 
to October 24th, 2006. 

For the photometric redshift study, we consider only fields 
with observations in all five photometric bands. Our parent sam- 
ple is therefore reduced to u*,g',r', i',z' images, namely 

- the 4x 1 deg^ Deep fields with longer integration times than 
the previous T0003 release (~ O.lmag deeper) used in 106, 
and 

- the 35 deg" of the Wide field that are completed in all filters 
and distributed as follows: 19 deg^ in Wl, 5 deg^ in W3 and 
1 1 deg2 in W4. 



' http : //terapix . iap . fr/cplt/oldSite/ 
Descart/suitimarycflitlswide .html 
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CFHTLS Wl 
CFHTLS Dl 
WDS deep 



323319-041200 



J23319-050800 



323319-050400 



323319-070000 



22929-0412C 



22929-050800 ! 22539-050800 



22929-050400 



22929-070000 



22539-060400 



39 38 



22150-041200 



22150-050800 



22150-050400 



36 

Ro 



21800-041200 



21800-050800 



21410-041200 



21800-060400 



21410-050800 



CFHTLS W3 
CFHTLS D3 
DFEP2 



41154 + 523831 



41201+514251 



0555+523831 



. 5955 + 523831 



140016+514231 



222054 + 011900 



CFHTLS W4 
WDS f22 



222C 34+002300 



???rl54-003ion 



21706+011900 



'21706+002300 



71706-003100 



21318+011900 



21115-003100 



20930+002300 



20930-003100 



336 335 



333 
Ro 



332 331 



CFHTLS 02 
zCOSMOS 



Fig. 1. Sky coverage of the CFHTLS Deep fields (D1,D2,D3,D4) and Wide fields (Wl, W3 and W4) overlapped with spectroscopic 
surveys: WDS F02, F22 and DEEP2. 



The sky coverage of the complete parent data set is is shown in 
figure [T] 

For each deep field, Terapix produces two stacks per filter 
One comprises the 25% best seeing images and the other con- 
sists of the 85% best seeing images. Since both types of deep 
stack have seeing better than one arc-second and we are primar- 
ily interested in very deep galaxy catalogues, we used the deepest 
85% to compute the photometric redshifts. 

Table [T ] sumrn arises the data used for this study and also in 
lllbert et alT(l2006l) . 



2.1 .3. Production of T0004 catalogues 

Full details of the processing and content of the T0004 release 
are described in the Terapix pagefl Further details on calibra- 
tion and stack and catalogue production can also be found also in 
iMcCracken et alj (fc008) and in the Mellier et al (2005) explana- 
tory document In what follows we summarise briefly how the 
input catalogues used for the photometric redshifts have been 
produced. 



- http : //terapix . iap . fr/rubrique . php?id_rubrique=241 
^ http : //terapix . iap . fr/cplt/oldSite/Descart/ 
NewterapixdocT0O02 . pdf 
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Table 1. Summary of the data used in this study and in 106. 

Ilbert et al (2006). CFHTLS release T0003 



Field 



Center Position 



Deep-Dl 02''25'"59', -04°29'40" 

Deep - D2 10''00'"28', +02° 12'30" 

Deep-D3 14''19'"27', +52°40'56" 

Deep-D4 22''15"'31% -17°43'56" 



Area (unmasked) 
0.77 deg" 
0.69 deg^ 
0.83 deg' 
0.82 deg^ 



80% completeness limit (;^g) Spectroscopic data 



25.1 
24.9 
25.7 
25.3 



VVDS Deep 
no 
no 
no 



This study. CFHTLS release T0004. 



Field 


Center Position 


Area (unmasked) 


80% completeness limit (i^g) 


Spectroscopic data 


Deep-Dl 


02''25"'59', -04°29'40" 


0.78 deg^ 


25.3 




VVDS Deep 


Deep - D2 


10''00'"28', +02°12'30" 


0.80 deg2 


25.1 




zCOSMOS 


Deep - D3 


14'' 19^27', +52°40'56" 


0.83 deg^ 


25.9 




DEEP2 


Deep - D4 


22''15'"31', -17°43'56" 


0.82 deg^ 


25.5 




no 


Wide - Wl 


02''18'"00', -07°00'00" 


19 (15.73) deg2 


> 24.3 




VVDS Deep 


Wide - W3 


14'47'"54', +54°30'31" 


5 (4.05) deg2 


> 24.4 




DEEP2 


Wide - W4 


22''13'"18', +0ri9'00" 


11 (8.87) deg2 


> 24.4 




VVDS F22 



Pre-processing of raw images (masking bad pixels, remov- 
ing the overscan, subtracting the dark and the bias, flat field- 
ing and illumination correction) is performed by the Elixir 
pipeline at CFHT (iMagnier & CuillandreH2004 ). The data are 
then transferred to Terapix to produce the stacked images 
and the final catalogues. Figure |2] shows the flowchart of the 
T0004 release production. Each image is first examined by 
the Terapix QualityFITS image quality control tool. During 
the QualityFITS step (hereafter QFITS-in) a set of quality 
assessements is produced, all individual input images are in- 
spected and evaluated and a weightmap image as well as an 
input catalogue are produced. This catalogue will be used later 
for the astrometric calibr ation. It is produced using SExtractor 
dBertin & ArnoutsI 19961) with appropriate settings for saturation 
levels to eliminate any spurious objects which could lead to an 
incorrect computation of the flux re-scaling during the astromet- 
ric matching process. A QFITS-in web page summarises the in- 
spection and is used as an ID-card of each image. All QFITS-ed 
images are then graded "A", "B" or "C", after a visual inspection 
of each web-page. 

With the QFITS-in information available, Terapix selected 
the T0004 image sample by applying the following criteria: 

- Terapix class: A or B grades (all images within the CFHTLS 
specifications; grade B are images that may show minor 
problems or with specifications very close to limits), 

- Exposure time higher than 60 s., 

- Seeing less than 1.3", except for m* (less than 1.4") 

- Airmass less than 1 .7 

- Skyprobe value lower than 2.0 magnitudes (security limit) 

Rejected images will no longer be considered. Terapix then 
uses the QFITS-in catalogues of the remaining images to derive 
the astrometric and photometric calibrations for the relea se. 

The astrometric solution is computed using SCAMF0 jBertinl 

120061) . SCAMP first examines all image headers and then split the 
exposures into a series of astrometric contexts. Each context sep- 
arates blocks of observing epochs where the instrument focal 
plane is in a fixed and (almost) stable configuration. The detec- 
tions and positions of astrometric sources on MegaCam images 
are derived by the cross-identification of sources of the QFITS- 
in catalogue with the 2MASS astrometric reference catalogue. 
For the T0004 release, the source matching exploration radius is 

'' http : // terapix . iap . f r/sof t/ scamp 



set to 2" for all Deep and Wide fields. A polynomial distortion 
model is then derived by minimising a weighted quadratic sum 
of differences in positions between the 2MASS and the QFITS- 
in matched sources, and, internally, between different QFITS-in 
catalogues with overlap regions. Using this information SCAMP 
can then compute external and internal errors. 

For the Wide data, the astrometric solution is performed only 
once for each Wide field by using all image catalogues simulta- 
neously, regardless the filter and the epoch. All images of a given 
Wide field are then calibrated globally and in a homogeneous 
way. For the 85% Deep data, the solution must be computed dif- 
ferently. The number of observing runs produces too many astro- 
metric contexts which results in a matrix that cannot be inverted 
with the current Terapix computing resources. The Deep field 
images are then divided into five sub-samples, one for each fil- 
ter. In order to improve the consistency and the robustness of the 
solutions found for each filter, a common set of extra images sur- 
rounding the field and shifted by about 30' with respect to Deep 
centre positions is added to each sub-sample. The consistency of 
each solution has been checked afterwards. 

In both cases, the Deep and the Wide calibrations worked 
well. For T0004, the rms internal error of Wide and Deep as- 
trometric solution is 0.017" and the mean rms external error is 
0.21" in both directions. After inspection, if acceptable, the as- 
trometric solution is then written in the .head file attached to 
each image. 

SCAMP is also used to derive a relative photometric calibra- 
tion by minimizing a weighted, quadratic sum of magnitude dif- 
ferences on the overlap regions between images. Images flagged 
as "photometric" by CFHT are used as anchor points and their 
CFHT magnitude zero point is written in the . ahead file. For 
the photometric calibration, SCAMP minimizes the quadratic sum 
of magnitudes using the overlap region between images. The 
flux of non-photometric images are then re-scaled accordingly. 
Typical re-scaling amplitudes in T0004 are +0.02 magnitude. 
However, for some highly non-photometric images it may reach 
±0.50 magnitudes. For these extreme outliers, the re-scaling it- 
self becomes more and more unreliable and the observed error 
on the rescaling value can be significantly higher than the typi- 
cal ±0.02 value. For this reason all images with rescaling greater 
than 0. 15 magnitudes were removed from the input list. 

After this post- SCAMP selection process images are divided 
into filter and tile positions and images are then resampled and 
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Global Scamp D85 
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Stack image 
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Chi2 image 
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Fig. 2. Details of the data processing performed by Terapix. 

co-added with SWariS dBertin et alj|2002l) . An ASCII polygon 
mask is produced at each tile position which can be visualised 
using the "DS9" image display software. Each Wide stack is 
only composed of images centred at the tile position and the 
overlapping pixels of nearest neighbour tiles are not used. For 
each stack a 1 degxl' strip masks the boudaries of fields in or- 
der to drop regions with highly elongated PSF and low signal- 
to-noise. 

T0004 stacks are produced by a median filter and a 
Lanczos-3 interpolation kernel. All stacks have 19354 X 19354 
pixels of 0.186" (i.e. exactly 1x1 deg^) and have a mag- 
nitude zero point set to 30. For the Wide survey, a stack 
is produced at each centre position listed in the Terapix 
web pag43. A SExtractor catalogue is then produced for 
each stack that is used to run QualityFITS and produce a 
QFITS-out "ID-card" of each stack. QFITS-out as well as 
all quality control data of all stacks are publicly available at 

' http : // terapix . iap . f r/sof t/ swarp 
^ http : //terapix . iap . fr/cplt/oldSite/Descart/ 
suitimarycfhtlswide . html 



http : //terapix . iap . £r/article . php?idjrticle=724. 
They have been used to carry out visual inspections and the 
validation of the T0004 release. 

If g'-, r'- and band data are available, Terapix au- 
tomatica lly produces a "ch i2 image" based on these three 
stacks ( Szalay et alJ Il999l) . Then SExtractor is run in 
dual-image mode on both the chi2 image and each stacked 
image (/, g' , i' , as well as u* and z' , if any). All cata- 
logues contain parameter values for all quantities listed in 
http : //terapix . iap . fr/article . php?idjrticle=628. 
In this study, we use the merged (m, g' , r' , i', z' ,) catalogues pro- 
duced by Terapix that include a limited number of parameters 
(only MAG_AUTO magnitudes, for example) plus the E{B - V) 
value at each source position derived from dust map images 
Schlegel et al. (1998). 

After removing the masked areas, the effective field- 
of-view is 0.78, 0.80, 0.83 and 0.82 deg^ for the Dl, 
D2, D3 and D4, respectively and 15.73, 4.05 and 8.87 
for the Wl, W3 and W4, respectively. Effective areas for 
each individual tile in the Wide fields are given in the 
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http : //terapix . iap . f r/cplt/table_syn_T0004 . html 
table. 

All catalogues include flux with uncertainties. However, in 
this work we follow 106 and we scaled the SExtractor flux 
error estimates by a factor of 1.5 to compensate for the slight 
noise correlation introduced by the image re-sampling during 
image stacking. 

Finally, a series of post-processing analysis is then carried 
out in order to make quahty assessments for each stack and, 
globally, for the whole release. The QFITS-in, SCAMP, QFITS- 
out output files are part of the quality assessment data. More spe- 
cific control files are also created using the merged (u*, g', r' , /', 
z') catalogues, like stellar colour-colour plots and comparisons 
between the T0004 stellar photometry and the Sloan Digital 
Sky Survey (SDSSfl This quality control data are available at 
http : //terapix . iap . f r/cplt/table_syn_T0(S(S4 . html. 

The comparison between CFHTLS T0004 and SDSS 
photometry is made using stars in common in the magni- 
tude range 17 < AB < 21.1. The results are displayed at 
http : //terapix . iap . £r/article . php?id^rticle=713. 
The mean offset between CFHTLS and SDSS, Acfhtls-sdss? 
varies from stack to stack in the range -0.03 < Acfhtls-sdss < 
0.03 magnitude. A calibration problem was noticed on M*-band 
images obtained during the period from March 2006 until 
October 24th, 2006 showing a mean offset of 0.21 mag. with 
respect to SDSS. Following the Terapix table all M*-band values 
obtained during this period were corrected for the offset derived 
from the SDSS. In the few CFHTLS fields with no overlap with 
SDSS, we applied a correction of 0.2 mag. 

2.1.4. Field-to-field photometric stability 

The Terapix pipeline uses common sources in overlap regions 
between each contiguous MegaCam field to derive a photometric 
magnitude zero-point correction after a field-to-field flux rescal- 
ing. It is computed by SCAMP to produce a Wide survey as uni- 
form as possible over all Wl, W2, W3 and W4 fields. The recipe 
relies on the CFHT magnitude zero points and photometric flags 
written in the FITS header and in the observing logs, respec- 
tively. It is possible that some series of marginally photometric 
images produce a poor flux rescaling solution. We minimised 
this risk by removing photometric outlier images from the T0004 
sample during the post- SCAMP selection process (See Section 
12.1.31 ). However, some fields may still show residuals, in partic- 
ular at the borders of the survey where the field-to-field rescal- 
ing cannot be done using the four MegaCam boundaries of each 
stack. The comparison between CFHTLS and SDSS photometry 
confirms the uniformity is good to within ±0.03 mag. over the 
whole fields, but there are a few fields that are off by as much as 
-0.10 mag. 

This problem can be partially overcome by Le Phare. The 
code uses spectroscopic redshifts to correct the magnitude off- 
sets (hereafter systematic offsets) in the CFHTLS wide fields 
that overlap with spectroscopic fields. The estimate of the field- 
to-field magnitude variations and how it can affect the photomet- 
ric redshift accuracy is discussed in a forthcoming section. 



2.2. Spectroscopic data 

We used spectroscopic redshifts (spectro-z's) from the VVDS 
survey for the Dl and Wl fields. The VVDS data were obtained 



with flie Visible Multi-Object Spectrograph (VIMOSH) installed 
at the ESO-VLT. The deep spectroscopic sample VVDS-0226- 
04 was selected in the magnitude range 17.5 ^ Iab ^ 24.0 and 
has a median redshift of ~ 0.67 ( Le Fevre et al.li2()05a ). Keeping 
only sources with a confidence level in the redshift measurement 
greater than 97% (class 3 and 4), our parent VVDS-0226-04 
spectroscopic sample is composed of 3880 galaxies. We matched 
3276 galaxies in the Dl field and 3356 galaxies in the Wl field 
with i'^g < 24. 

We also used t he public spectroscopic redshifts from the 
zCOSMOS survey dLillv et al. 2007) which overlaps with the 
D2 field. The zCOSMOS-bright spectra were obtained with the 
VIMOS spectrograph and were selected at Iab ^ 22.5. 3915 
spectro-z's with a confidence level greater than 99% are usable 
for the D2 field. 

In addition, our spectroscopic calibration sample includes 
5936 spectro-z's fr om the third data relea se of the DEEP2 survey 
jDavis et al. l l2003l and lD"avis et al ] |2007l) . This sample overlaps 
with the D3 field and covers a small area in the W3 field (see 
FiglU- Spectra were obtained by the spectrograph DEIMOSQ 
mounted on the Keck II telescope. As for the VVDS sample we 
used the most secure spectro-z's, with quality flag 3 or 4, cor- 
responding to confidence level greater than 95%. The DEEP2 
data taken in the Extended Groth Strip region have been prese- 
lected in the range 18.5 < Rab ^ 24.1, with a selection based 
on color and surface brightness only aimed at maximizing the 
number of galaxies over the number of stars (in contrast with the 
rest of the DEEP2 survey where selection criteria to target only 
higher-redshift galaxies are used). A smaller subsample of 310 
spectro-z's are in the W3 field. 

Finally, the VVDS- Wide survey (VVDS-F22 field, 
iGariUi et alj (|2008|) ) overlaps with the W4 field. The VVDS- 
Wide covers a 4 deg^ contiguous area in the W4 field. The 
spectro-z catalogue includes 11 228 galaxies, 167 type I AGNs, 
and 6 748 stars. Using the most secure spectro-z (confidence 
level greater than 97%, class 3 or 4) our final sample comprises 
3 854 galaxies in the magnitude range 17.5 < Iab ^ 22.5. 

D4 is the only field without any spectroscopic coverage. 



3. Photometric redshifts for the CFHTLS-T0004 
catalogue 

3.1. Photometric redshift computation with Le Phare 

Photometric redshifts were computed using a standard 
template-fitting procedure and calibrated using spectroscopic 
redshifts. We use d t he code Le PharA^ described in 
Arnouts e t all d 1999b and lArnouts et alj ( |2002|) with the addition 
of the optimisation procedure presented in 106. 

We first selected a set of reference SED templates. These 
were the same CFHTLS -optimised templates as 106 used. 

The original set of templates is composed of four observed 



galaxy spectra from lColeman et aTj(ll98Q) (hereafter C W W) and 
two starburst galaxy spectra from iKinnev et al.l (Il996h . The ul- 
traviolet {A < 2000A) parts of the spectra have been linearly 
extrapolated and the near infrare d parts are extrapolated us ing 
the synthetic models proposed bv iBruzual & CharlotI (l2003l) . It 
is clear that the large variety of galaxy spectra observed in the 



http : //www . sdss . org/data 



http : //www . eso . org/ sci/ f acilities/paranal 
/instruments/vimos/overview.html 
' http : //loen. ucolick . org//Deimos/deimos .html 
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Universe cannot be represented by only a small number of op- 
timised templates. However, adding too many templates creates 
degeneracies between observed colors and redshifts. For this rea- 
son we use only a small number of spectral types. 

The six templates were "optimised" for the CFHTLS us- 
ing 2867 spectroscopic redshifts from the VVDS deep survey. 
The optimisation procedure, described in 106, consists of "blue- 
shifting" the observed photometric data to rest-frame using the 
VVDS spectroscopic redshifts. A low resolution optimised ver- 
sion of the template is obtained by averaging the rest-frame 
fluxes of all the galaxies (as shown in Fig. 5 of 106). Finally, the 
templates were linearly interpolated between spectral types in 
order to cover the full redshift-color space. In total, 62 optimised 
templates are generated. 

Each SED template was then redshifted onto a grid of in- 
terval 6z - 0.04 and convolved with the filter transmission 
curves (including instrument efficiency). The photometric red- 
shifts were derived by determining which SED template pro- 
vides the best match to the observed colours (minimisation of the 
merit function). The galaxy internal reddening, E(B - V), was 
included as a free parameter in the template-fitting procedure. 
The values allowed for E{B - V) were derived from the Small 
Magellanic Cloud extinction law dPrevot et alj 119841) . varying 
from to 0.2 for Scd and later types, and no reddening was al- 
lowed for earlier types. 

Since no near infrared data were available to break the 
colour-redshift degeneracies between z < 0.5 and z > 2.5 (in 
general caused by an inability to distinguish between between 
the Balmer and Lyman breaks), a "prior" on the redshift distribu- 
ti on has be en applied following the Bayesian approach described 
in Benftez' (^2000). We used the parameterised redshift distribu- 
tion estimated from the VVDS deep spectroscopic survey. In 106, 
the redshift distribution was estimated by spectral type and as a 
function of apparent magni tude . On the b asis of the luminosity 
functions measured by Ilbe rt et all (l2005h we reject unrealistic 
objects with Mg> < -24. 

A redshift Probability Distribution Function (PDFz) is com- 
puted for each object using the merit function, PDFz oc 
exp{-x^{z)/2). The PDFz is measured every 6z = 0.04. The best 
redshift is estimated via a parabolic interpolation of the PDF. If 
a second peak is found in the PDF with a height larger than 5% 
of the main peak, the corresponding redshift is given as a second 
solution. 

In addition to the best derived from the galaxy library 
(hereaft er X^ai )' ^ best x^ computed using the stellar library of 
iPicklesI (Il998h is derived for each object (hereafter ;t'^,^,.). Both 
X^'s are use to compute the star-galaxy classification as will be 
explained in a later Section. 

3.2. Correction of systematic offsets 

As dem onstrated by iBrodwin et"aLl ( |2006l ), 106 and lllbert et al.l 
( |2008|) . systematic off'sets are often found between the best-fit 
SED templates and the observed apparent magnitudes in a given 
filter. Uncertainties in the zero-point calibration of the photomet- 
ric data as well as imperfect knowledge of galaxy SEDs or of 
filter transmission curves are responsible for these offsets. They 
may produce additional biases in the photometric redshift mea- 
surements and must be corrected. 

In this work, we recomputed these offsets to account for the 
changes in our spectroscopic and photometric catalogues with 
respect to those used in 106. These changes are: 



- The T0004 catalogues have been regenerated after a com- 
plete re-processing of all data;. 

- The CFHTLS Wide data were not used by 106. These new 
stacks cover 35 deg^ and combine many MegaCam tiles; 

- New spectroscopic samples are now available for many 
CFHTLS fields (which were not available in 106). It is now 
possible to compute offsets for many new fields and to derive 
field-to-field systematic offsets. 

Following the procedure used in 106, we minimised the SED 
template ;(f^'s at a fixed spectroscopic redshift with an additional 
free off'set per band that account for the systematics. The offsets 
are then derived iteratively. For all spectroscopic surveys, the 
calibration was performed using the most secure spectroscopic 
redshifts. Stars were not used for the offset computation. 

In this study we limited the spectroscopic samples to /^^ < 
22.5. As shown by 106 this correction depends only weakly on 
the spectroscopic limiting magnitude. Each offset correction was 
derived and applied independently for each field covered with 
spectroscopic data (see section IZ21 ). Table |2] shows the system- 
atic offsets derived for the Deep and Wide fields in each band. 
The systematic offsets vary from -0.041 {g') to 0.077 (m*) mag- 
nitudes and show a small dispersion between the fields, on the 
order of 0.01 mag. For this reason we added 0.01 magnitudes 
in quadrature to the error estimate to account for the systematic 
offset uncertainties. 

D4 is the only field which does not overlap with any spec- 
troscopic survey. For this field, we used a mean offset correc- 
tion, computed using a combined catalogue of Dl/VVDS deep, 
D2/ZCOSMOS and D3/DEEP2. 

3.3. Pliotometric/spectroscopic comparison for tine T0004 
Deep fields 

We evaluate in the two following sections the photometric 
redshift accuracy by comparing photometric redshifts against 
spectro-z's. We first consider the Deep fields; in the following 
Section we consider the Wide fields. 

As in 106, we define two relevant quantities: 

- the photometric redshift dispersion, crA2/(i+;j. It is defined 
using the normalised median absolute deviation o-az/{1+z,) - 
l.48xmedian{\Az\/{l+Zs)), a robust approximate of the stan- 
dard deviation, 

- and the outlier rate, rj, also called catastrophic errors, defined 
as the proportion of objects with |Az| > 0.15 x (1 + Zj), 

where Zp is the photometric redshift, Zs the spectroscopic red- 
shift, and Az = Zp - Zs- 

For the VVDS Deep and DEEP2 surveys we can make com- 
parisons to i'^g < 24. In the D2 field, for the brighter zCOSMOS 
survey comparisons cannot be made fainter than /^^ - 22.5. 
Figure[3]shows the comparison between photometric and VVDS 
spectroscopic redshifts for the Dl field; zCOSMOS for the D2 
field and DEEP2 for the D3 field. The dispersion (0-4^/(1+2 )) for 
Dl and D3 in the range 17.5 < i\g < 24 is 0.028 and 0.030 re- 
spectively, and the dispersion for D2 in the range 17.5 < i\g < 
22.5 is 0.026. The outlier rate (?/) is 3.57% in Dl, 3.62% in D3 
and 1 .35% in D2. These results demonstrate that the photometric 
redshift accuracy is comparable for the three deep fields. Slightly 
better results are found in D2 as this sample has a brighter limit- 
ing magnitude. We checked this assumption by cutting the lim- 
iting magnitudes for the Dl and D3 at 17.5 < /^^ < 22.5. We 
found the same dispersion for the three fields, crA„/(i+;,) - 0.026, 
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Table 2. Systematic offsets of i'^g < 22.5 limited samples for the Deep and Wide fields. 
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Fig. 3. Photometric redshift accuracy in the Deep survey. We compare photometric redshifts with spectro-z's for D 1 with WDS deep 
(left), for D2 with zCOSMOS (middle) and for D3 against DEEP2 (right). The brighter zCOSMOS survey is limited to i'^g = 22.5. 
In each panel we show the dispersion, cr a;/(i+,j = 1.48 x median{\Az\/{l +Zs)), and the outlier rate, rj, the proportion of objects with 
|Az| > 0.15 X (1 H- Zj). In the top panel, triangles represent objects for which we found a second solution in the PDF. In the bottom 
panel, the object density in the zphot-zspec plane is plotted. 

as well as slightly different but better outlier rates for the three ing the photometric redshift accuracy or increasing the fraction 
fields: 1.68forDl, 1.35 for D2 and 2.36 for D3. of outliers. 



The consistent results found for these three deep fields are 
reassuring. Although our templates have been optimised using 
one single field (Dl) with one set of spectra (WDS Deep), they 
can be used with confidence on other deep fields without degrad- 
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3.4. Photometric/spectroscopic comparison for the T0004 
Wide fields 

The photometric redshift accuracy is function of the signal- 
to-noise ratio (hereafter S /N) of the photometric data (e.g. 
Margoniner and Wittman 2007). Since the Wide survey is much 
shallower than the Deep survey the median S /N is lower at a 
given magnitude in the Wide than in the Deep (see Figure |4] for 
the Dl and Wl fields). Therefore, we expect less accurate photo- 
metric redshifts for the wide in comparison with the deep fields 
at the same magnitude. 

First, we measured the photometric redshift accuracy as a 
function of decreasing S /N (i.e. increasing magnitude) in the 
wide fields. 

Noise is estimated using magnitude errors from 
SExtractor's MAG_AUTO, A/w, rescaled by a factor 1.5 
(See Section 2.1.3). Using the VVDS deep spectroscopic 
sample, we measured the photometric redshift precision to 
i'^g = 24. In Table [3] we show the dispersion and the outlier rate 
for the Wl field as a function of limiting magnitude. We found 
that cr A7/(i+7,) is larger at fainter magnitudes. At i'^g < 24 the 
dispersion is 0.053. However the number of outliers increases 
dramatically beyond i'^^ - 23, reaching 10% at i'^^ < 24. 
We conclude that photometric redshifts are reliable down to 
i'j^g = 22.5 in the Wide field (with a median S /N - 33). Fainter 
than i'^g = 23 the number of outliers becomes important. 

Next we compared the photometric redshift accuracy be- 
tween the different Wide fields. We limited the comparison at 
i'^g < 22.5 (median S /N = 33) which is the limiting magni- 
tude of the zCOSMOS and VVDS-Wide spectrocopic surveys. 
Figure |5] shows the comparison between photometric and spec- 
troscopic redshifts (Wl against VVDS deep, W3 against DEEP2 
and W4 against VVDS Wide). The dispersion crA;/(i+7^) is found 
to be identical for all fields, ranging from 0.036 to 0.039 (0.037 
for Wl, 0.036 for W3 and 0.039 for W4). The outlier rate, //, 
is 2.81%, 3.55% and 3.79% for Wl, W3 and W4, respectively. 
As for the Deep survey, these results show that one can success- 
fully use templates optimised in an independent field to reach 
comparable accuracy. 

We also noticed that despite being calibrated with a small 
spectroscopic sample, the photometric redshift accuracy and out- 
lier rate in the D3 field are comparable to the other fields. For 
the purpose of photometric redshift calibration, it appears that 
for the same amount of time dedicated to a spectroscopic survey 
it is more important to cover a larger area with bright spectro- 
scopic redshifts (which can correct for inhomogeneities in the 
photometric calibration) rather than a smaller. 

Finally, in order to perform a fair comparison between the 
Deep and Wide photometric redshift one must consider sam- 
ples with equal S IN. Therefore, in order to match the median 
S /N - 43 corresponding to i'^g < 24 in the Deep, we limit the 
magnitude to i'^g = 22 in the Wide (from figure |4]i and measure 
a dispersion of 0.032, similar to the value of 0.028 found for the 
Dl photometric redshifts. 

3.5. Effects of tlie tile-to-tile systematic offset variatiori 

From FiglUone can see that CFHTLS fields have only partial 
spectroscopic coverage. The offset calibration can only be per- 
formed in sub-areas where there is spectroscopic coverage, and 
we assume the systematic offsets measured in a sub-area of the 
field can be applied to the whole field. 

This assumption is most valid for the CFHTLS deep fields 
since each spectroscopic sample fits within one tile, but more 



80 






60 






40 






20 




Wide ^ — -^^^ 





Deep — — - 



21.5 22.0 22.5 23.0 23.5 24.0 24.5 25.0 

iob < 

Fig. 4. Median signal-to-noise ratio with respect to the limiting 
magnitude in the Wide (Wl) and the Deep (Dl) field. The S /N 
is 1.0857 /Am, taken in the /' band. The dashed line represents 
the median S /N - 33, corresponding to the value of the Wide 
sample cut at /^g < 22.5. 

Table 3. Photometric redshift accuracy in the Wl and Dl fields 
for several limiting magnitudes. 
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questionable for the CFHTLS wide fields which are ten times 
larger and are composed of many continuous tiles. As it was dis- 
cussed in Section 2.1.4, the Aotcfhtls-sdss photometric com- 
parison performed by Terapix shows the differences with SDSS 
have a mean scatter of +0.03 magnitudes which corresponds to 
the photometric scatter from tile to tile. Table |2] confirms these 
trends in our calibrations. 

Only 15 of 35 tiles in the Wide fields have spectroscopic data 
(not all centered in the tile). Therefore, the calibration of the sys- 
tematic offsets cannot be performed tile by tile in the wide fields. 
Instead, they are averaged over several tiles (Wl has four over- 
lapping tiles over 19 with VVDS deep, W3 has two overlapping 
tiles in five with DEEP2 and W4 has nine overlapping tiles in 1 1 
with VVDS wide). Therefore, tile-to-tile photometric calibration 
variations are an additional source of error which degrades the 
photometric redshift accuracy in the Wide fields. 

We assessed the impact on the photometric redshift accuracy 
of this additional source of error, assuming a tile-to-tile variation 
of 0.03 magnitudes (see Sect. 2.1.3). We simulated variations of 
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Fig. 5. Photometric redshift accuracy in the Wide survey at i'^^ < 22.5. We compared photometric redshifts with spectro-z's for Wl 
with WDS deep (left), for W3 with DEEP2 (middle) and for W4 with WDS Wide (right). In each panel we show the dispersion, 
Caz/(1+z,) - 1-48 X median(\Az\/(l +Zj)), and the outlier rate, 77, the proportion of objects with |Az| > 0.15 x (1 +Zs)- In the top panel, 
triangles represent objects for which we found a second solution in the PDF. Bottom: density of objects in the zphot-zspec plane is 
plotted. 



the photometric calibration on fields with spectroscopic data as 
follows: 



We repeated this procedure 1000 times and computed the mean 
and standard deviation for the photometric redshift dispersion cr, 
and for the outfier rate, 77. 



1 . we computed the systematic offsets from a field with spec- 
troscopic data (with unchanged photometry); 

2. we modelled the tile-to-tile zero-point variation with a 
Gaussian distribution with zero mean and a dispersion of 
0.03. We applied offsets to this field, but on the 5 bands inde- 
pendently (i.e. we drew an independent value for each band); 

3. we computed the photometric redshifts but we kept the offset 
corrections derived at the first step. Since we did not compute 
the systematic offsets again this sequence is equivalent to 
derive systematic offsets from one field and apply them on 
another field; 

4. finally, we estimated the uncertainties (photometric redshift 
dispersion and outlier rate). 



Figure |6] shows the results for the Deep and Wide fields. As 
expected, photometric redshift accuracy degrades. On average, 
o"Az/(i+:,) is 21% worse in the Deep field and 14% worse in the 
Wide field, whereas they were selected to have comparable S IN. 
The worse degradation observed for the Deep field is likely due 
to the numerous faint objects, more sensible to photometric vari- 
ations. The outlier rate shows a similar trend, with a value 13% 
worse in the Deep sample and 17% in the Wide sample. This 
analysis shows the need for spectroscopic surveys spread over 
large areas. 
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Fig. 6. Degradation of the photometric redshift accuracy if the 
offset correction is not computed on the same tile. The results 
show the variation of the outlier rate (left panels) and photomet- 
ric redshift accuracy (right panels) if the photometry varies of the 
order of 0.03 magnitude. Top panels: Dl Deep field at i'^g < 24 
(compared with 0-^2/(1+;,) - 0.028 and 77 = 3.57%, represented 
by the red dashed lines, with no variation added in the photome- 
try). Lower panels: Wl Wide field at i'^g < 22.5 (compared with 
Caz/ci+z,) — 0.037 and rj - 2.81%, represented by the red dashed 
lines, with no photometric perturbation applied). 
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Fig. 7. Photometric redshifts for the Wide field (Wl, y-axis) 
compared with photometric redshifts for the Deep field (D 1 , x- 
axis). The sample is selected at 17.5 < /^^ < 22.5. We fit in 
the range [0.2,0.7] to better measure remaining systematic dif- 
ferences. These are small compared to the dispersion and we 
conclude there is a good consistency between the two samples. 



3.6. Wide/Deep comparison 

As described in section IZTl Wide and Deep fields have been re- 
duced independently. The systematic offsets were also calibrated 
independently. For this reason, galaxies common to a Deep and 
a Wide field may have two different photometric redshift values. 

To check if results from Wide and Deep surveys are consis- 
tent, we performed a photometric redshift/photometric redshift 
comparison using overlapping sources in the Wl and Dl fields. 
We selected 14484 common sources in Dl and Wl in the range 
17.5 < i'j^g < 22.5. 

The figure Q shows the Wl photometric redshift versus the 
Dl photometric redshift. We find the dispersion between the 
wide and deep photometric redshift samples is o-a;/(i+z,) - 
0.026, whilst the outlier rate is 77 = 2.71 similar to values found 
for the Dl sample, demonstrating that the Wide and Deep pho- 
tometric redshifts are well consistent with each other 



4. Star/galaxy selection 

Bright samples {i'^g < 22.5) can be highly contaminated by stars. 
For example in the CFHTLS W4 field stellar sources comprise 
up to 50% of sources and can be proble matic for analy s es sen - 
sitive to such mis-classifications. As in iSchultheis et al.l (l2006h . 
we used the VVDS spectroscopic data to improve the star/galaxy 
selection. 

Thanks to purely flux limited selection of the VVDS spec- 
troscopy survey, we can test our ability to separate stars and 
galaxies at the depth of the Wide and Deep CFHTLS surveys. 
We have derived an efficient technique to separate stars from 



galaxies based on both morphological and multi-colour infor- 
mation. 

4. 1 . Size selection 

A widespread technique to separate stars from galaxies involves 
comparing the size of sources with the local shape of the PSF. 
Bright objects comparable in size to the local PSF are identi- 
fied as stars while those with larger size are selected as galax- 
ies. Figure |8] illustrates the difference in size between galaxies 
and stars for two tiles of the Wl and W4 fields. The size of the 
sources is characterised by the half-light radius, r/,, which is de- 
fined as the radius which encloses 50% of the object flux. The 
lower panels of figure [8] show the size (r/,)-magnitude (/^^) di- 
agram (iFahlman et alJ ll994). Due to statistical measurement er- 
rors of r/, and to the large PSF variation over the MegaCam field 
of view the dispersion of the r/, distribution is broad. 

Spectroscopic redshifts can be used to securely distinguish 
stars and galaxies and to clearly define star-galaxy selection cri- 
teria calibrated on bright sources. We define a "true" star as an 
object with z, = and a "true" galaxy an object with a secure 
spectroscopic redshift Zs > 0. Figure [8] shows the positions of 
spectroscopically confirmed stars (asterisks) and spectroscopi- 
cally confirmed galaxies (black squares) in the r/, vs /^^ diagram 
for two tiles in the wide fields. At bright magnitudes, a strin- 
gent cut in r/, can reliably separate stars from galaxies. We then 
set this cut on each CFHTLS Deep field and Wide tile indepen- 
dently: at i'^g < 21 the r/, distribution is predominantly com- 
posed of stars. 
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The histogram of r/, values is close to a normal distribution, 
slightly skewed towards larger r/, by the galaxies (top panel of 
Figure |8]i. For this reason we fit the r/, distributions at /^^ < 21 
of all Deep and Wide tiles independently by a Gaussian. We de- 
note and cr^j the best fit mean and standard deviation values 
respectively. In the Gaussian approximation of the r/, distribu- 
tion, almost all stars should lie in the range < r/, < ju,,^ + 3cr,.^ - 

^hlimit' 

4.2. Colour selection 

A purely size-based selection as described in the previous sec- 
tion limits reliable star-galaxy classification to only the brightest 
sources. As shown in Figure |8] the proportion of spectroscop- 
ically identified galaxies fainter than i'^g -21 which satisfies 
the size criterion r/, < r/,;,„„-, increases towards fainter magni- 
tudes, where galaxy and star populations become increasingly 
mixed. Adding an additional colour-based selection can improve 
the star-galaxy separation at fainter limits. 

Following the galaxy spectral type classification based on the 
galaxy template fitting criteria (namely, j^'^^,^^.^,), we also charac- 
terise each source by its stellar spectral type, based on a fitting 
of st ellar ternplate s (namely, xlmr)- We used a set of templates 
from lPicklesI ( Il998i) to estimate xltar- 

Ideally, a galaxy would satisfy the relation x\ai ^ xltm- 
Yioviewer, x'g^ig^^. and;^^^^ were derived from independent SED 
libraries and with a different number of parameters. The rel- 
evance and the statistical significance of the comparison be- 
tween xlaiaxy ^"'^ X^tar ^sve been assessed by using the spec- 
troscopic sample. We used a free parameter to account for the 
differences between each estimate and found that the crite- 
rion < xltarl'^ is robust and can be applied to all sources up 
to = 23. 

When only optical data {u*,g', r', i',z') are used, stellar and 
galaxy colors overlap (shown in figure |9] is the distributions 
for stars and galaxies in the spectroscopic sample). As in the 
case of a purely size-based selection, a selection based solely on 
a multi-colour criterion cannot be used in isolation. The most ro- 
bust star-galaxy separation consists in using both the;^'^ estimate 
and the size criterion. Each source is therefore defined by three 
parameters: its size, r/,, its galaxy index, X^aiaxv' stellar 

indexxiar 

4.3. Quality of the star/galaxy classification 

In order to assess the accuracy of the method we tested our se- 
lection against a strict spectroscopic selection. We defined the 
incompleteness of the parent galaxy sample as the number of 
galaxies lost after the selection compared to the number of galax- 
ies and the star contamination of the final galaxy sample as the 
number of stars misidentified as galaxies compared to the num- 
ber of estimated galaxies. 

At i'^g < 21, stars are purely selected with their sizes. At 
fainter magnitude, we defined a galaxy as an object with r/, > 
fir^ + 2>a-ri, orx^ai < xliai-l'^ ("or" is used to have a galaxy sample 
as complete as possible). Stars are defined as r/, < fir,^ + 30"^^ and 
xltar < -^^gfl/- W^ stopped the selection at i'^g > 23, where all 
objects are flagged as galaxies. 

Table|4]shows the results in the Dl, Wl and W4 fields. The 
star contamination is strongly reduced even for the fields where 
the proportion of stars is more than 50%. 
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Fig. 9. x^ distributions for "true" stars and "true" galaxies (ac- 
cording to the VVDS Wide spectroscopy) in the field W4. Top: 
all objects. Middle: objects larger with r/, > r/,;,,,,,-,. Bottom: ob- 
jects with r/i < rhumit- if used alone, the x^ selection or the r/, 
selection lead to misidentification. On optimal way to separate 
stars from galaxies is to combine the two pieces of information. 

The methods works the best with a small "PSF", leading to 
a small r/,;,„„Y (< 2.5 pixels). 

The colour-based selection is helpful to keep the galaxy 
sample as complete as possible where the "PSF" is large, e.g. 
rhiimii ~ 3 pixels (the incompleteness being 10% instead of 
3%, otherwise). 

4.4. Conclusion 

The need for an accurate star/galaxy selection has arisen from 
the W4 field, where numerous stars contaminate the galaxy sam- 
ple. We showed that the star/galaxy separation is simple for 
bright objects (i'^g < 21) but becomes less reliable when faint, 
small objects are mixed together We used the spectroscopic 
sample to build up and to assess a reliable estimator for the na- 
ture of an object, using both the size and color of objects. We 
succeeded in keeping an incompleteness always below 3% (in 
the worse "PSF" case) while we significantly reduced the star 
contamination (from 50% down to 7% in average). In the best 
cases (small "PSF" and a few stars in the field), the method 
reaches 1 . 1 % of star contamination with less than 1 % of incom- 
pleteness. 

5. Photometric redshift analysis 

From this Section onwards we consider only "reliable" photo- 
metric redshifts. These are defined as redshift of a source esti- 
mated from five photometric bands located in unmasked regions 
and which fulfilled passed the size-SED star-galaxy classifica- 
tion described in the previous Section. The source must be in the 
following magnitude range: 17.5 < i'^g < 24 for the Deep fields 
and 17.5 < i'^g < 22.5 for Wide fields. Hereafter we focus on the 
reliable photometric redshift samples to assess the quality red- 
shift catalogues, to derive error estimates and compute the the 
redshift distribution in both the Wide and Deep fields. 
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Fig, 8. Star/galaxy separation in one tile of the Wl field, 022539-041200 (left panel), and one from the W4 field, 221706+002300 
(right panel). Top: r/, distribution and the corresponding gaussian fitting. The red line is the + 3cr^,^ cut. Bottom: i'^g vs r/, in 
pixels. Squares and asterisks are VVDS-Wide spectroscopically identified stars and galaxies respectively. Both fields have different 
specifications; Wl overlaps with a deeper spectroscopic sample and has a low stellar density, whereas W4 has a larger proportion 
of stars but a brighter spectroscopic limit. 



Table 4. Star/galaxy selection results, the table shows the galaxy (PDF) of each object. By determining the redshifts where 

incompleteness (Inc.) and the star contamination (Cont.) for the Ax^(z) - X^(z) ~ Xmi„ - 1' we determine a low (ziefi) and a 

Wl field, with a large PSF size and the highly star crowded field, high (Zright) redshift value corresponding to a 68% confidence 

W4. "Stars" represents the intrinsic proportion of stars in the interval. 

field. We define the error estimated by "Le Phare" cr,^ as follows: 



Deep < 24) 

Field /■;,;,„„, Stars Inc. Cont. 
■^1 277 nWc 080% Y7[% 



Wl (i'^g < 22.5) 



Field 


hiimil 


Stars 


Inc. 


Cont. 


022539-050800 


3.1 


23.5% 


3.1% 


0.4% 


022539-041200 


2.8 


18.8% 


2.3% 


0.9% 


022929-041200 


3.1 


20.6% 


2.9% 


0.7% 


AH 


3.0 


19.0% 


2.6% 


1.0% 


W4 


d',, < 22.5) 






Field 


hiimit 


Stars 


Inc. 


Cont. 


221318-003100 


2.2 


57.2% 


0.6% 


7.7% 


221318+002300 


2.1 


52.5% 


0.0% 


9.0% 


221318+011900 


2.6 


44.7% 


0.8% 


2.5% 


221706-003100 


2.7 


53.2% 


1.0% 


7.9% 


221706+002300 


2.3 


51.9% 


1.0% 


6.4% 


221706+011900 


1.9 


51.7% 


0.4% 


7.3% 


222054-003100 


2.5 


52.6% 


1.0% 


6.1% 


222054+002300 


24 


52.8% 


1.2% 


7.3% 


222054+011900 


2.1 


45.8% 


0.3% 


6.3% 


All 


2.3 


51.8% 


0.8% 


6.9% 



5.1. Error estimates 

"Le Phare" computes an estimate of the photometric redshift 
uncertainty from the redshift Probability Distribution Function 



|zfe/,(68%)-z„v»(68%)| 

o"- = ■ (1) 

~" 2 

The reliability of this error estimate is examined as follows. 
We first compare cr,^ to the variance of the difference between 
spectroscopic and photometric redshifts AzWe then use cr,^ to 
assess the accuracy of the photometric redshifts over the entire 
range of magnitudes and redshifts. 

Figure [To] shows the distribution of differences between the 
spectroscopic and the photometric redshifts (Az), for the wide 
Wl, W3 and W4 fields (left) and for the deep Dl and D3 (right) 
fields. On this Figure is overlaid a Gaussian with a mean of zero 
and a standard deviation equal to the dispersion derived from the 
data, that is o-A;._wiV/f/(i+-,) - 0.038 and o-Aj;_D(,f/)/(i+z,) = 0.029. 
A very good agreement is found in the Deep and Wide fields 
demonstrating that the errors are nearly Gaussian. The distri- 
bution shows an extended tail and departs from the Gaussian 
approximation for large errors which is caused by the outliers. 
The distribution shows a slight bias in the Wide field only of 
Azsys = median{{zp - z.s)/(l + Zs)) - -0.005. A closer look at 
the bias as a function of the magnitude shows sUghtly underesti- 
mated photometric redshifts for very bright objects (/^^ < 20.0) 
and overestimated photometric redshifts at fainter magnitude 
(see table |5]l, but this bias is always less than 1%. 

The cumulative distributions of the estimated errors com- 
pared to the photometric/spectroscopic error can also be com- 
puted. Assuming a Gaussian errors distribution we expect that 
68% of the objects have a photometric redshift within the range 
Zs ± lo"jp, recalling that cr,^_ is the uncertainty estimated by 
"Le Phare". Figure [TT] shows the cumulative distribution of 
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Fig. 10. Histogram of Az for the Wide (W1,W3 and W4, left) and the Deep (D1,D3, right) fields. A Gaussian function with a mean 
of and a dispersion corresponding to the photometric redshift dispersion is superimposed. In both cases the error distribution is very 
well represented by this Gaussian function. 



Table 5. Median photometric redshift bias in magnitude slices 
in the field Wl. 



T03 - T04 error comparison 







17.5-20.0 


0.0070 


20.0-20.5 


-0.0016 


20.5-21.0 


-0.0055 


21.0-21.5 


-0.0080 


21.5-22.0 


-0.0078 


22.0-22.5 


-0.0073 



kp - ■Z.!l/(lo"zp- From the two panels, we find that 71.3% of 
sources in the Wide and 70.8% in the Deep satisfy the relation 
kp - Zs\ < icr^^, showing the reliability of the Gaussian approxi- 
mation of error estimates. For the remainder of this work photo- 
metric redshift errors will be derived from cr, of "Le Phare". 

-/J 

We can use this error model to investigate the accuracy of 
photometric redshifts over the complete magnitude and redshift 
ranges of our catalogues. Figure [12] shows the fraction of photo- 
metric redshifts in the Deep survey with cr-._ < 0.15 x (1 -H Zp) 
as function of i'^^ apparent magnitude. For photometric redshifts 
without outliers one would expect to find almost all objects in 
this range. As seen from the comparisons with spectroscopic 
redshifts a small number of objects have catastrophic errors at 
z < 0.2 and z > 1.5 where degeneracies in the colour-template 
space becomes important. 

The solid and dotted lines show the photometric redshift er- 
rors for the galaxies in common between T0003 (106) and T0004 
(this work). Both studies show comparable results with a slight 
improvement in T0004 for z > 0.2. 

The comparison between photometric redshift errors in Wide 
and Deep fields is shown in FigurjT3] 95% of sources have an 
error within +0.15 x (1 -H Zp), for i'^g < 24 in the Deep field and 
i'^g < 22.5 in the Wide fields. This result is consistent with the 
comparisons with spectroscopic redshifts where the outlier rate 
never exceeds 5%. 
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Fig. 12. Photometric redshifts errors as function of magnitude 
and redshift for the Deep fields for this study (solid line) and for 
106 (dotted line). 



The dependence of photometric redshift errors on redshift is 
illustrated in figure[T4l As already pointed out in 106, the photo- 
metric redshift errors in the Deep survey is lowest for z < 1.5. 
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Fig. 11. Cumulative errors for the Wide (W1,W3 and W4) and the Deep (D1,D3) fields. As expected for a Gaussian error distribution, 
the photometric redshifts fall within zs ± Icr^^ 68% of the time where Icr, , is the estimated error 



5.2. Redshift distribution 

The redshift distributions of the T0004 Deep and Wide sources 
were derived from the histogram of photometric redshifts. 

We took into account three sources of uncertainties in the 
redshift distribution: the uncertainties on the photo-z, the cosmic 
variance and the Poisson errors. 

We derived the photo-z uncertai nties on the redshift bins as 
follows. Following IFu etal.1 (|2008|) . we model the PDF by a 
normalised Gaussian distribution with nns cr-^ - |z/f/f(68%) - 
Zng/ir(68%)|/2 and mean Zp for each galaxy. From each individ- 
ual PDF we randomly drew a redshift and built the histogram of 
the redshift distribution for the whole sample. We repeated the 
process 100 times and computed the dispersion for each redshift 
bin. 

Cosmic variance and Poisson errors are combined into one 
single error value, where the field-to-field variance is estimated 
directly from the data. In the Deep field, we simply computed the 
field-to-field scatter over the four independent fields and divided 
it by Vi. 

In the Wide field, the effective area is different for Wl, W3 
and W4, so we estimated the field-to-field variance in a different 
way. We cut Wl into compact subareas and computed the field- 
to-field variance as function of their angular size (see Fig. ]re- 
fcv). We then extrapolated the results to derive the field-to-field 
variance corresponding to the full size of each Wide field. Due 
to a correlation between adjacent subareas, the field-to-field vari- 
ance is likely understimated. As for the Deep field, we divided 
by V3 the error, assuming the 3 Wide fields are independent. 

Finally we fit the redshift distribution by an ana- 
lytic function. We use the parameterised form proposed by 
[VanWaerbeke et al. (2001): 

niz)^^—l^]\J-l^J] (2) 
ZQr[^)\zol [ \zol ) 

and we fit over the range 0.0 < z < 1 .5 for the Wide redshift 
distribution and over the range 0.0 < z < 2.5 for the Deep one. 



Table 6. Mean and median redshifts in the Deep and Wide fields, 
for several limiting magnitudes. Both values are computed di- 
rectly from the data. 





Deep 






Wide 




'ab < 


< z > 


Zmedian 


'as < 


< z > 


Zmedian 


20.5 


0.36 


0.34 


20.5 


0.37 


0.35 


21.0 


0.40 


0.37 


21.0 


0.41 


0.39 


21.5 


0.46 


0.42 


21.5 


0.46 


0.44 


22.0 


0.51 


0.47 


22.0 


0.51 


0.49 


22.5 


0.58 


0.54 


22.5 


0.56 


0.55 


23.0 


0.65 


0.62 


23.0 


0.62 


0.60 


23.5 


0.74 


0.70 


23.5 


0.68 


0.65 


24.0 


0.84 


0.77 


24.0 


0.73 


0.70 


24.5 


0.93 


0.83 


24.5 


0.77 


0.75 


25.0 


1.00 


0.88 


25.0 


0.81 


0.70 



Figure [16] shows the redshift distribution derived for the 
fields Wl, W3, W4 and Dl, D2, D3, D4 grouped in a Wide and 
Deep redshift distributions, respectively. As expected, the cos- 
mic variance is lower in the Wide field which covers 35 deg^. 

We estimated the mean redshift of the Deep and Wide fields 
for several limiting magnitude. Table|6]shows the mean and me- 
dian redshift computed directly from the data for several lim- 
iting magnitudes. As expected, the mean and median redshifts 
increase as the limiting magnitude of the sample increases. 

The top panel of Figure [17] compares the Wide and Deep 
redshift distributions and their fit for a same magnitude selection 
at 17.5 < i'^g < 22.5. The redshift distributions of the Deep and 
Wide fields agree well, when a same magnitude limit is adopted. 

Conclusions 

We computed photometric redshifts for the release T0004 of the 
Deep and Wide CFHT Legacy Survey produced by Terapix. In 
this work, only T0004 fields observed in all filters (m*, g' , r', 
/', and z') relevant for photo-z studies have been selected. Our 
photometric catalogues cover four independent Deep fields of 
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Fig. 14. Estimated photometric redshift error for the Wide (Wl) 
and Deep (D 1 ) fields as function of redshift. Errors are binned in 
redshift with an interval = 0.04, which explained the steps in the 
distribution. Few objects only are in the range z > 2 in the Wide 
field so errors are less reliable in this range. 
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Fig. 13. Photometric redshift errors as function of magnitude and 
redshift for both the Wide and the Deep fields (upper panel). 
Lower panel: A zoom of the redshift range 0.2 < z < 1 .5. 



1 deg^ each (Dl to D4) and 35 deg^ distributed over three inde- 
pendent Wide fields (W1,W3 and W4). 

We define "reliable" photometric redshifts as redshifts de- 
rived for sources classified as galaxies according to a joined size- 
SED star-galaxy classification criterion, located in unmasked re- 
gions with five-band photometric data and with / < 1000. 



These objects have a magnitude in the range 17.5 < i'^^ < 22.5 
and 17.5 < i'^g < 24 for the Wide and Deep samples, respec- 
tively. 

In total, the photometric redshift catalogues contain 244,701 
reliable redshifts in the Deep field within the magnitude range 
17.5 < < 24 and 592,891 reliable redshifts in the Wide fields 
within the range 17.5 < /^^ < 22.5. 

Our photometric redshift method is based on a SED template 
fitting procedure using the code "Le_Phare". Following Ilbert et 
al (2006) spectroscopic redshifts were used to optimise our SED 



Fig. 15. Relative field-to-field variance of subareas in Wl as 
function of the area (overplotted is the Poisson error). Results are 
extrapolated to derive an estimate of the total field-to-field vari- 
ance in the Wide field. Uncertainties on the field-to-field vari- 
ance were derived from the data using a Jackknife estimator 



templates and to correct the systematic offsets between the ob- 
servations and SEDs. A Bayesian prior is used to break the color- 
redshift degeneracies. 

A total of 16,983 spectroscopic redshifts from the VVDS 
Deep, VVDS Wide, DEEP2 and zCOSMOS spectroscopic sur- 
veys provide a spectroscopic coverage for each field (except D4). 
The large quantity of spectroscopic data enabled an independent 
systematic offset correction for each field. A mean correction 
was applied to D4. These offsets vary only a small amount be- 
tween the fields, on the order of 0.01 mag, demonstrating the 
excellent field-to-field stability of the CFHTLS photometry. 

We estimated the accuracy of our photometric redshifts by 
comparing them to spectroscopic catalogues. The Deep field 
photometric redshifts were compared to the VVDS and DEEP2 
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Fig. 16. Redshift distribution in the Deep (top) and Wide (bot- 
tom) fields. Error bars comprise the photo-z error and the field- 
to-field variance (cosmic variance and Poisson error). Error bars 
show the uncertainty for the full Deep field (top),computed from 
the field-to-field variance between D1,D2,D3 and D4 and multi- 
plied by a factor 1 / V4, since all 4 Deep field are independent. 
In the Wide field, the field-to-field variance were computed for 
several subareas and extrapolated to the size of the Wide fields 
and finally divided by V3, as the three Wide fields are also inde- 
pendent. 

spectroscopic samples limited at 17.5 < /^^ < 24. We found 
a stable dispersion of 0.028-0.030 and an outlier rate of about 
3-4%. For the Wide fields, the dispersion is 0.036-0.039 and 
the outlier rate is about 3-4% in the range 17.5 < /^^ < 22.5 
(comparing with VVDS Deep, VVDS Wide, zCOSMOS and 
DEEP2). The systematic bias between the input spectroscopic 
samples and the photometric redshifts derived with Le Phare 
keeps below 1% over 17.5 < i'j^g < 22.5 for the Wide fields and 
no significant bias is found in the Deep fields. We noticed that 
the dispersion is constant or slighlty increasing for S /N > 40, 
with a more significant increase at lower S /N values. The out- 
lier rate increases dramatically in the Wide beyond i'^^ = 23 (5% 
and 10% at /^^ < 23 and < 24). Our results are comparable 
to similar analyses done in the same CFHTLS Wide fields by 
[Erben et al. (2008) and Brimioulle et al. (2008). 

Each CFHTLS Wide field is composed of several 1 deg^ 
tiles. Unfortunately, the spectroscopic coverage is not sufficient 
to perform a calibration of the systematic offsets for each tile. 
We investigated the effect of tile-to-tile photometric offsets on 
the accuracy of our photometric redshifts. Taking an uncertainty 
of 0.03 magnitudes for the photometric calibration in each band, 
we found that the photometric redshifts can be degraded up to 
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z 

Fig. 17. Distribution of the photometric redshifts (i'^g < 22.5) 
computed on the new data release (T0004) for the Wide fields, 
plotted separately on the lower panel and plotted combined on 
the upper panel. 

14% in the Wide and 21% in the Deep. Future wide-field spec- 
troscopic surveys such as VIPERS will allow to reduce this prob- 
lem. 

Some catalogues in the Wide fields contain 50% of stars at 
i'^^ < 22.5. A reliable galaxy/star separation is important to sci- 
entifically exploit these catalogues. We combined both the size 
and multi-colour data to perform a robust star/galaxy separation. 
The star selection is purely based a size critea for the brightest 
sources (i'^g < 21). Since some galaxies are unresolved at fainter 
magnitudes, we added multi-colour information to avoid reject- 
ing galaxies. Using the spectroscopic data to assess the quality 
of our classification, we found that our criteria reduces the stel- 
lar contamination from 50% (in W4, the field with the highest 
stellar density) to 7%, while keeping a galaxy sample more than 
98% complete. Star contamination is as low as 1-2% in some 
other fields, like Wl or Dl. 

Finally, we used the spectroscopic redshifts to evaluate the 
68% error estimate computed by "Le Phare" for the photomet- 
ric redshifts. The error distribution, found to be very well ap- 
proximated by a Gaussian, was in excellent agreement with 
the real photo-z dispersion. Confident with the 68% error esti- 
mate for each photometric redshift, we determined the accuracy 
of our photometric redshifts over a larger magnitude and red- 
shift. Our photometric redshifts are most accurate in the range 
0.2 <z< 1.5. 

The redshift distributions for the Deep and Wide fields have 
been modeled using empirical formula of IVan Waerbeke et al.l 
(l2001h . The distributions and the Zmedian and mean < z > de- 
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rived from the Deep and Wide samples are consistent with 106, 
Fu et al.' ('2008') as well as the recent CFHTLS-Wide studies of 
Erben et al. (2008) and Brimioulle et al. (2008). 

The photometric redshifts calculated here will be an es- 
sential tool in realising the full potential of the CFHTLS sur- 
vey. In particular, cosmological parameter estimation, galaxy- 
galaxy lensing, galaxy biasing studies and investigation of the 
halo occupation function for galaxies will greatly benefit from 
this homogeneous and well-calibrated set of photometric red- 
shifts. Additionally, upcoming large spectroscopic surveys like 
VIPERS will benefit from the increased observing efficiency that 
pre-selection using photometric redshifts can offer 

Our CFHTLS T0004 Deep and Wide photometric redshifts 
are made available to the community through the TerapixOH and 
CENCOSQ databases. 

Acknowledgements 

We thank the CFHT, Terapix and CADC for their assistance 
and considerable work in producing the CFHTLS data. We also 
acknowledge the VVDS consortium for providing the spectro- 
scopic redshift catalogues. This work has made use of the VizieR 
catalogue access tool, CDS, Strasbourg. We thank S. Colombi, 
T. Erben, H. Hildebrandt, B. Menard, S. Seitz and L. van 
Waerbeke for useful discussions. MK is supported by the CNRS 
ANR "ECOSSTAT", contract number ANR-BLAN-0283-04. 
We acknowledge the CNRS-Institut National des Sciences 
de rUnivers (INSU) and the French Progamme National de 
Cosmologie (PNC) for their support for the CFHTLS. This re- 
search was supported by ANR grant "ANR-07-BLAN-0228". 



Ilbert, O., Capak, P., Salvato, M., et al. 2008, Ai'Xiv 0809.2101 

Ilbert, O., Tresse, L., Zucca, E., et al. 2005, A&A, 439, 863 

Kilbinger, M., Benabed, K., Guy, J., et al. 2008, ArXiv 0810.5129 

Kinney, A. L., Calzetti, D., Bohlin, R. C, et al. 1996, ApJ, 467, 38 

Le Fevre, O., Guzzo, L., Meneux, B., et al. 2005a, A&A, 439, 877 

Le Fevre, O., Vettolani, G., Garilli, B., et al. 2005b, A&A, 439, 845 

Lilly. S. J., Le Fevre, O., Renzini, A., et al. 2007, ApJS, 172, 70 

Magnier, E. A. & Cuillandre, J.-C. 2004, PASP, 1 16, 449 

McCracken, H. J., Ilbert, O., Mellier, Y., et al. 2008, A&A, 479, 321 

Mobasher, B., Capak, P, Scoville, N. Z., et al. 2007, ApJS, 172, 117 

Parker, L. C., Hoekstra, H., Hudson, M. J., van Waerbeke, L., & Mellier, Y. 2007, 

ApJ, 669,21 
Pickles, A. J. 1998, PASP, 110, 863 

Prevot, M. L., Lequeux, J., Prevot, L., Maurice, E., & Rocca-Volmerange, B. 

1984, A&A, 132, 389 
Rowan-Robinson, M., Babbedge, T., Oliver, S., et al. 2008, MNRAS, 386, 697 
Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 525 
Schultheis, M., Robin, A. C., Reyle, C, et al. 2006, A&A, 447, 185 
Szalay, A. S., Connolly, A. J., & Szokoly, G. P 1999, AJ, 117, 68 
Tereno, I., Schimd, C, Uzan, J.-P, et al. 2008, ArXiv 0810.0555 
Van Waerbeke, L., Mellier, Y, Radovich, M., et al. 2001, A&A, 374, 757 
Wolf, C, Meisenheimer, K., Rix, H.-W., et al. 2003, A&A, 401, 73 



References 

Amouts, S., Cristiani, S., Moscardini. L., et al. 1999, MNRAS, 310, 540 
Arnouts, S., Moscardini, L., Vanzella, E., et al. 2002, MNRAS, 329, 355 
Ball, N. M., Branner, R. J., Myers, A. D., et al. 2008, ApJ, 683, 12 
Bem'tez, N. 2000, ApJ, 536, 571 

Benjamin, J., Heymans, C, Semboloni, E., et al. 2007, MNRAS, 381, 702 

Berge, J., Pacaud, F, Refregier, A., et al. 2008, MNRAS, 385, 695 

Bertin, E. 2006, in Astronomical Society of the Pacific Conference Series, Vol. 

351, Astronomical Data Analysis Software and Systems XV, ed. C. Gabriel, 

C. Arviset, D. Ponz, & S. Emique, 112-+ 
Bertin, E. & Amouts, S. 1996. A&AS, 117, 393 

Bertin, E., Mellier, Y, Radovich, M., et al. 2002, in Astronomical Society of the 
Pacific Conference Series, Vol. 281, Astronomical Data Analysis Software 
and Systems XI, ed. D. A. Bohlender, D. Durand, & T. H. Handley. 228-+ 

Boulade, O., Chariot, X., Abbon, P., et al. 2000, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 4008, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. 
M. lye & A. F Moorwood, 657-668 

Brimioulle, F, Lerchster, M., Seitz, S., Bender, R., & Snigula, J. 2008, ArXiv 
0811.3211vl 

Brodwin, M., Brown, M. J. I., Ashby. M. L. N., et al. 2006, ApJ, 651, 791 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Coleman, G. D., Wu, C.-C, & Weedman, D. W. 1980, ApJS, 43, 393 

Davis, M., Faber, S. M., Newman. J., et al. 2003, in Presented at the Society 

of Photo-Optical Instiumentation Engineers (SPIE) Conference, Vol. 4834, 

Discoveries and Research Prospects from 6- to 10-Meter-Class Telescopes II. 

Edited by Guhathakurta, Puragra. Proceedings of the SPIE, Volume 4834, pp. 

161-172 (2003)., ed. P Guhathakurta, 161-172 
Davis, M., Guhathakurta, P, Konidaris, N. P, et al. 2007, ApJ, 660, LI 
Erben, T., Hildebrandt, H., Lerchster, M., et al. 2008, ArXiv 081 1.2239 
Fahlman, G., Kaiser, N., Squires, G., & Woods, D. 1994, ApJ, 437, 56 
Fu, L., Semboloni, E., Hoekstra, H., et al. 2008, A&A, 479, 9 
Garilh, B., Le Fevre, O., Guzzo, L., et al. 2008, A&A, 486, 683 
Hildebrandt, H., Wolf C, & Benitez, N. 2008, A&A, 480, 703 
Ilbert, O., Amouts, S., McCracken, H. J., et al. 2006, A&A, 457, 841 



" http://terapix.iap.fr 

'2 http://cencos.oamp.fr/CFHTLS 



